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ABSTRACT 



With high-resolution infrared data becoming available that can probe the formation of high- 
er ' mass stellar clusters for the first time, appropriate models that make testable predictions of 
these objects are necessary. We utilize a three-dimensional radiative transfer code, including a 
hierarchically clumped dusty envelope, to study the earliest stages of super star cluster evolution. 
We explore a range of parameter space in geometric sequences that mimic the hypothesized 
evolution of an embedded super star cluster. The inclusion of a hierarchically clumped medium 
can make the envelope porous, in accordance with previous models and supporting observational 
evidence. The infrared luminosity inferred from observations can differ by a factor of two from 
the true value in the clumpiest envelopes depending on the viewing angle. The infrared spectral 
energy distribution (SED) also varies with viewing angle for clumpy envelopes, creating a range 
in possible observable infrared colors and magnitudes, silicate feature depths and dust continua. 
General observable features of cluster evolution differ between envelopes that are relatively opaque 
or transparent to mid-infrared photons. For optically thick envelopes, evolution is marked by 
a gradual decline of the 9.8/im silicate absorption feature depth and a corresponding increase 
in the visual/ultraviolet flux. For the optically thin envelopes, clusters typically begin with a 
strong hot dust component and silicates in emission, and these features gradually fade until the 
mid-infrared PAH features are predominant. Epr the models with a smooth dust distribution, 
the Spitzer MIPS or Herschel PACS [70]- [160] color is a good probe of the stellar mass relative 
to the total mass, or star formation efhciency. Likewise, the IRAC/MIPS [3.6]-[24] color can be 
used to constrain the Ri„ and Ron* values of the envelope. However, dumpiness confuses the 
general trends seen in the smooth dust distribution models, making it harder to determine a 



1. INTRODUCTION 

Super star clusters (SSCs) are massive young 
star clusters with high stellar densities (> 10'^ 
stars pc^'^) that form under extreme pressures, 
often found in merging galaxy systems, galactic 
nuclei, and blue compact dwarf galaxies (Whit- 
more 2002, and references therein). With masses 
typically in excess of lO^M©, they are the most 
massive type of stellar cluster known. Studies of 
the most nearby analogues Westerlund 1 (Clark 
et al. 2005) and R136 (Bosch et al. 2009) show 
that these objects are consistent with expectations 
that they are globular cluster progenitors; this 
expectation is strengthened by numerous stud- 
ies of SSC populations outside the Local Group 
(e.g. Schweizer et al. 1996; Holtzman et al. 1992; 
O'ConneU et al. 1994; Conti & Vacca 1994; John- 
son et al. 1999). Because of their high density of 
massive stars, models show that SSCs have the 
ability to disperse metals from supernovae ejecta 
to large distances, trigger further star formation 
episodes, and act as a launching mechanism for 
super-galactic winds (Wiinsch et al. 2008; Murray 
et al. 2010). Therefore SSCs can have a major 
impact on their host galaxies. 

Tracking the early evolution of SSCs requires 
long wavelength observations due to the dust- 
enshrouding birth envelopes that surround them 
at young ages. Embedded SSCs were first detected 
in the radio (e.g. Kobulnicky & Johnson 1999; 
Turner et al. 2000) in low-metallicity blue com- 
pact dwarf galaxies (BCDs) as compact sources 
with high ionizing fluxes. Numerous other exam- 
ples have also been found in the infrared, in star- 
bursts such as Arp 220 (Shioya et al. 2001) and 
again in BCDs (Sauvage & Plante 2003). 

Radio observations have confirmed the masses 
of embedded SSCs to be in the range found for 
their optical counterparts. Furthermore, the cal- 
culated high ionizing fluxes are equivalent to hun- 
dreds or thousands of 0-stars and electron densi- 
ties are ^ 10^~^cm~^ (Beck et al. 2002; Johnson 
et al. 2003; Johnson & Kobulnicky 2003). These 
densities are similar to those found in Galactic HII 
regions, but arc significantly lower than typical 
densities associated with Active Galactic Nucleus 
Broad Line Regions (> lO^cm"^) suggesting that 
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SSCs are also heated by young stars in a starburst 
environment, and not AGN. 

Despite advances made this last decade in un- 
derstanding embedded SSCs as starburst events, 
their early evolution is still poorly understood. 
Properties such as the mass and physical size of 
the clouds from which they are formed and the 
amount of gas that is turned into stars (i.e. the 
star formation efficiency; Ashman & Zepf 2001) 
are not well-constrained. In addition, the most re- 
cent radio, near-infrared, and ultraviolet data on 
these objects suggest that the thick envelopes are 
porous, allowing a significant fraction of UV light 
to leak from the system (Reines et al. 2008; John- 
son et al. 2009; Thuan et al. 2005). This project 
is a first attempt at modeling the infrared spec- 
tral energy distributions of embedded super star 
clusters. By investigating a large variation in the 
values of the input parameters in a hierarchically 
clumpy and porous envelope, these models can be 
used to constrain the dusty envelope geometry and 
star formation efficiency of embedded super star 
clusters. 

The decision to model embedded SSCs in the 
infrared is based on the fact that they are most 
visible in the infrared and radio; at these young 
ages, hot stars heat the dust in the circum-cluster 
envelope and excite free-free radio emission inside 
the Hil region. Therefore these wavelengths are 
essential for studying the dust properties during 
the embedded phase, the dust mass, and the enve- 
lope geometry. The radiative transfer models pre- 
sented in this work offer predictions about what 
we expect embedded SSCs to look like in the in- 
frared and what wavelength observations offer the 
best diagnostic capability. With telescopes like 
JWST on the horizon, its high spatial resolution 
and mid-infrared sensitivity making it far better 
than existing infrared space facilities (Sonneborn 
2008), we may finally be entering a period when 
qiicstions about initial dust mass and geometry of 
a SSCs embedded phase can be answered on a 
global scale. In the near term, observations using 
Herschel and Spitzer may be used to probe unre- 
solved embedded SSCs, and these models can be 
used to estimate their physical size and geometry. 

In this paper we present simulated infrared im- 
ages, SEDs, and colors of embedded SSCs along 
a geometric sequence that mimics the evolution 
of a young embedded super star cluster. We dis- 
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cuss the models in detail in §2, then present model 
images, SEDs and colors in §3. We compare un- 
resolved populations to the models in §4, discuss 
the limitations of these models when comparing to 
resolved populations in §5, and conclude in §6. 

2. MODELS 

2.1. The Radiative Transfer Method 

The models presented in this paper calculate 
radiative transfer from thermal dust grains and 
stochastic small grains, and scattering from the 
dust. The code simulates a three-dimensional ge- 
ometry, and is based on previous work on thermal 
dust grains in two dimensions (see Whitney et al. 
2003; Chakrabarti & Whitney 2009), updated to 
include three dimensions and clumpy dust distri- 
butions (Indebetouw et al. 2006), and stochastic 
emission from polycyclic aromatic hydrocarbons 
and very small grains (PAHs/VSGs; Wood et al. 
2008). 

The code uses a Monte-Carlo technique that is- 
sues 'photon packets' into a dusty envelope from 
a central source. Packets that encounter large 
grains either change the temperature in the grid 
space in which the encounter occurs, according to 
the prescription of Lucy (1999), or are scattered 
according to a modified Henyey-Greenstein func- 
tion (Cornette & Shanks 1992). A small num- 
ber of packets encounter PAHs/VSGs, and are 
re-emitted according to emissivity templates from 
Draine & Li (2007); see Wood et al. (2008) for 
details about emission from PAHs/VSGs in our 
models. 

Embedded SSCs are believed to have porous 
envelopes that allow a significant portion of the 
ultraviolet (UV) light to escape from the system 
(supporting observational evidence is provided in 
§2.5 from Thuan et al. 2005; Reines et al. 2008; 
Johnson et al. 2009). To account for this fact, 
we have included a hierarchically clumped density 
structure for the dusty envelope, using the pre- 
scription of Elmegreen (1997). 

In order to confine the parameter space inves- 
tigated in this model set, we have varied only four 
input parameters that most influence the output 
spectral energy distribution: the inner and outer 
dust radii in a spherical geometry; the mass in 
stars divided by the total cluster mass (i.e. the 
'star formation efliciency' or SFE; Ashman & Zepf 



2001); and the fraction of dust that is smoothly, 
as opposed to fractally, distributed. The many 
other parameters that remain flxed, such as the 
central star cluster mass, the IMF, the central 
source luminosity and age, the dust species used, 
the volume fractal dimension adopted, and the ra- 
dial dust distribution exponent, are all described 
in §2.2 to §2.4. 

2.2. The Dust Composition 

We use the standard dust model size distribu- 
tion derived in Kim et al. (1994) that fits an ex- 
tinction curve with Ry = 4 using the prescription 
of Cardelli et al. (1989). The size distribution is 
not a simple power law, but is derived from a Max- 
imum Entropy Method solution, which gives the 
best solution compatible with the extinction data. 
The chosen Ry value has been found to be ap- 
propriate for the more dense regions of molecu- 
lar clouds (Whittet et al. 2001) and therefore we 
make the assumption that it will also fit embedded 
SSCs. 

We adopt the dielectric functions of astronomi- 
cal silicate and graphite grains from Laor & Draine 
(1993). Our grain model includes a layer of water 
ice on the grains, covering the outer 5% in ra- 
dius. The ice mantle increases the opacity of the 
grains at all wavelengths, but the change is most 
pronounced longward of 35/im. The slope of the 
opacity with wavelength is the same in the infrared 
for grains with and without the ice mantle, so any 
error introduced by the inclusion of an ice man- 
tle is a scalar factor at the wavelengths of concern 
in this paper. The ice dielectric function is from 
Warren (1984). 

The PAHs/VSGs are templates from Draine 
& Li (2007). The energy density of the radia- 
tion field relative to the interstellar radiation field 
(ISRF) can vary between 0.5 < f/ < 10^ and the 
ionization fraction is that used in Draine & Li 
(2007) to constrain the models to reproduce the 
observed Milky Way spectrum. Only grains that 
are < 200A in size are considered small grains. 
Grains larger than this are assumed to thermalize 
after photon interactions, and are therefore con- 
sidered large grains. The PAH/VSG mass fraction 
in dust used in this paper is 3.55%, which is the 
average value found for 65 galaxies in the SINGS 
survey (Draine et al. 2007). This is therefore an 
acceptable value to expect for extragalactic em- 
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bedded SSCs. 

2.3. The Central Source 

In these models, the central cluster is treated 

as a point source with the spectrum of a 10® Mq 
star cluster, produced by the Starburst99 spectral 
synthesis models (Leitherer et al. 1999). The input 
metallicity is solar and the initial mass function is 
a Salpeter IMF (Salpeter 1955); changing these 
parameters was found to have little effect on the 
overall shape of the resulting infrared SED, vali- 
dating the decision to hold them constant. The 
age is set to 1 Myr throughout the geometric se- 
quence to ensure that the UV luminosity remains 
constant. Since evolution is expected to take place 
in less than three or four million years, prescribing 
a fixed age is warranted - see §2.5 for a discussion 
about the evolutionary pace for embedded super 
star clusters. The instantaneous starburst mode 
was chosen, without continuum emission. The lu- 
minosity of the source is 1.6 x lO^L©. 

2.4. The Envelope Properties 

The dust and gas around the central source 
is distributed spherically with variable inner and 
outer radii and envelope mass. The initial en- 
velope mass is set for star formation efficiencies 
(SFEs) of 5%, 10%, 15%, 25%, and 50%, which 
is a large range of possible values based on ob- 
servations of star-forming regions (Bonnell et al. 
2010; Murray et al. 2010). The dust can be dis- 
tributed fractally and smoothly, and the percent- 
age of dust that is smoothly distributed is varied 
from between 100% smooth to 1% smooth (i.e. 0% 
to 99% clumpy). 

The models mimic a cluster's evolution by mov- 
ing the inner radius out towards the outer radius 
in a geometric sequence; Table 1 lists the sequence 
for the three possible outer radii we have chosen, 
and Figure 1 shows how the geometric progression 
roughly follows what is expected for a real SSC. 
At very young ages, each star will realistically in- 
habit its own compact Hii region, but as the clus- 
ter evolves, these compact regions will grow and 
merge into a single, large cavity. Assuming a sta- 
tionary stellar wind as in the semi-analytical and 
numerical models presented by Tenorio-Tagle et 
al. (2007), this cavity will grow and disperse, even- 
tually leaving the central SSC exposed. The ini- 



tial inner radius of the models, correlated to intra- 
cluster dust, is O.lpc, and moves out to the outer 
radius, which remains fixed. 

The SFE values from 5% to 50% presume ini- 
tial cocoon masses in the range 1.9 x IO^Mq down 
to 1 X IO^Mq. For a given SFE value and outer 
radius, the average mass density is constant with 
radius, so that the mass of the shell in a geomet- 
ric sequence scales with the inner and outer radii 
as: M oc R^^f — i?f„. The average mass densities 
range from 10^ cm~^ at the densest to 10 cm~^ 
at the least dense, where the densest models are 
for a SFE of 5% and an outer radius of 25pc while 
the least dense models are for a SFE of 50% and 
an outer radius of lOOpc. 

The envelope is assumed to have a standard 
dust-to-gas ratio of 0.01. It is furthermore seeded 
with PAHs/VSGs to a depth equal to Ay = 1.5 
from the ionizing source. This is because the 
UV flux is sufficiently attenuated after this depth 
that we expect very little UV heating of small 
grains in the rest of the envelope. Comparisons 
between models run with PAHs/VSGs through- 
out the dusty envelopes and models PAHs/VSGs 
in to a depth oi Ay = 1.5 show that PAH/VSG 
emission changes only slightly for models where 
Ay < 250. 

There are five clumpy dust fractions presented: 
0, 0.1, 0.5, 0.9, & 0.99 clumpy. The spread in 
clumpy fraction is the complete range of conceiv- 
able clumpy values that could be found. The vol- 
ume fractal dimension was chosen as the best value 
from Elmegreen (1997) which matched the turbu- 
lent intercloud medium, D=2.3. The fractal length 
L is taken to be 3.792 so that the maximum den- 
sity of the clumps can be computed per hierarchic 
level, p{h) = L(^-D)h+HD ^^ere h is the hierar- 
chical level in question and H is the maximum 
hierarchical level, which we have set equal to five; 
higher H values create fractal sizes smaller than 
the resolution elements in our three-dimensional 
grid. 

2.5. Supporting Observational Evidence 
for the Model Parameter Choices 

The model parameters were chosen to be con- 
sistent with observations of embedded SSCs, and 
the parameter space is corroborated with observa- 
tional evidence from the literature in this section. 
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Fig. 1. A cartoon illustrating the early stages of a super 
star cluster's evolution and our corresponding model esti- 
mates. The smallest inner radius presented in this paper, 
Rin = 0.1 pc, is likely to correspond to intra-cluster dust in 
a real cluster (left-most figure) , which will only exist while 
the constituent stars still have individual cocoons. 

With regard to dust grains types, the PAH 
models from Drainc & Li (2007) arc the most suit- 
able available to date. Ongoing research (Draine 
& Li 2007; Galliano et al. 2008) that accounts for 
different PAH ionization fractions and redder ex- 
citing sources than those used here would likely 
change individual PAH band fluxes. However, 
since we are not interested in PAH ratios in this 
work, the inclusion of these processes does not im- 
pact the results of this study. 

The ice mantle on the dust grains is suitable for 
dense, cold regions, such as young star formation 
regions, where a large part of the dust is not being 
heated by the embedded stars and the water-ice 
cannot be sublimated from the grains. At any op- 
tical depth it has been shown that the water-ice 
features will be masked and therefore hard to de- 
tect (Robinson & Maldoni 2010). Since there is 
observational evidence for water associated with 
compact, young star clusters (Brogan et al. 2010) 
and the slope of the opacity functions with and 
without water ice are the same, we have left the 
dust with ice mantles throughout the entire geo- 
metric sequence. 

Fixing the cluster age is based on upper limits 
of the time it takes for a cluster to emerge from its 
envelope. Observations with HST show optically 
visible clusters with ages of just a few Myr (e.g. 
Johnson et al. 2000). Radio observations have also 
suggested the same approximate age, so we can 
assume a fiducial emergence timescale of about 
three Myr for a super star cluster (Whitmore 2002; 
Kobulnicky & Johnson 1999). Since the UV con- 



tinuum, which most affects the resulting infrared 
SED, does not significantly change between ages 
of and 3-4 Myr, setting the cluster age to 1 Myr 
throughout the sequence is warranted in order to 
limit the number of free parameters. 

The models are run under the assumption that 
the central star cluster is a point source and the 
envelope is perfectly spherical; this is a necessary 
simplification to confine the parameter space in 
this study. However, the high stellar densities 
and small half-light radii discovered for optically- 
visible SSCs (Tacconi-Garman et al. 1996; de 
Marchi et al. 1997) suggest that the central source 
is fairly localized in its formation envelope. Fur- 
thermore, placing dust to within O.lpc of the cen- 
tral point source will produce the high dust tem- 
peratures expected from intracluster dust in a real 
embedded SSC. Dust sublimation is not a concern 
even for these most embedded model phases: the 
dust sublimation radius, assuming a sublimation 
temperature of 1600K, is O.OSpc, well within the 
smallest inner radius value adopted. 

The outer radii of 25, 50, and lOOpc were cho- 
sen based on estimates from observational results. 
Vacca, Johnson, & Conti (2002) model infrared 
and radio observations of 'ultra-dense' Hii regions 
in Henize 2-10 as scaled-up versions of Galactic 
ultra-compact Hll regions to derive their radial 
size. Hunt et al. (2005) based their estimates on 
models of global dwarf galaxy SEDs for which the 
radii of the actively star-forming regions are de- 
rived based on infrared luminosities and tempera- 
tures. 

Lastly, in order to be consistent with the grow- 
ing body of evidence that suggests that the ISM 
is clumpy in a large variety of environments, a 
clumpy envelope was included around the cen- 
tral star cluster. That the ISM is not uniform is 
well established (e.g. Verschuur 1995; Elmegreen & 
Efremov 1997; Kim et al. 2008). In star formation 
environments, recent work has begun to show that 
the same is true. Thuan et al. (2005) shows that 
the molecular hydrogen in a BCD galaxy is clumpy 
because emission is visible in the near-IR, which 
is sensitive to dense H2, but not in the UV, which 
is sensitive to diffuse H2. Reines et al. (2008) and 
Johnson et al. (2009) show that for different clus- 
ters, more than 40% and 50% of the UV flux is de- 
tected outside the embedded sources. This could 
either be because the ISRF is very strong around 
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the clusters or because there is light leaking from 
the clusters. Most recently, model-dependent re- 
sults by Eldridge & Relaiio (2010) shows that 
about 50% of the ionizing flux is leaving the gi- 
ant Hii region NGC604, though this latter result 
is highly dependent on ages derived from observa- 
tions of Wolf-Rayet stars. 

3. RESULTS 

3.1. Simulated Images 

Images of the clusters at several stages along 
the geometric sequence are shown in Figure 2 in 
near-IR (J, H, and K bands), mid-IR {Spitzer 
IRAC 3.6/im, 4.5/^m, and 8.0//m bands; Fazio ct 
el. 2004), and far-IR {Spitzer MIPS 24^m, TO^im, 
and 160/im bands; Rieke et al. 2004) light from left 
to right. They were created from the SFE — 10% 
models with an outer radius of 50pc, clumpy dust 
fraction of 0.9 and inner radii of 1, 5, 20, and 45pc 
from top to bottom. 

One of the striking features of these simulations 
is that a significant amount of the near-IR light 
can escape the clumpy envelope, as seen by the 
scattered light on the inside of the cloud surfaces. 
Due to the scattered light, the near-IR images ap- 
pear blue in color, because scattering favors bluer 
photons (in this case, J-band photons are scat- 
tered more regularly than K-band photons). The 
mid-IR band images, sensitive to warm dust and 
PAH/VSG emission, appear red because of the 
dominant PAH feature and high dust continuum in 
the 8/im band. As the inner radius moves outward, 
the MIR colors become more red because the inner 
dust temperature becomes lower, removing what 
little hot dust radiation in the small-wavelength 
bands there was. Additionally the hole appears 
to be the brightest feature in the mid-IR because 
the shell is being lit up by the central source. The 
far-IR bands probe both sides of the infrared dust 
peak, with the 24/im emission showing hotter dust 
and 70 and 160/im bands showing cold dust emis- 
sion. Since the dust peak is typically near 70/im, 
these images therefore appear green. 

3.2. The Change in Ultraviolet Flux with 
Clumpy Dust Fraction 

Motivated by observations suggesting that 
roughly 50% of the UV flux leaks from embed- 
ded star clusters (see Reines et al. 2008; Johnson 




Fig. 2. — Simulated color images for a geometric (pseudo- 
evolutionary) sequence of a natal SSC for a SFE of 10%, 
outer radius = 50 pc, and clumpy fraction = 90%. Images 
are shown for an inner radius of the cocoon of 1 pc, 5 pc, 
20 pc, and 45 pc, from top to bottom respectively. The col- 
ors shown are for combinations of J, H,and K (left), IRAC 
3.6/jm, 4.5/im, and 8fim (middle), and MIPS 24/im, TO/tm, 
and 160/im (right). 

et al. 2009, for details on the observations), we 
plot the fraction of UV light lost as a function of 
the Rin/Rout ratio for the 99% clumpy models in 
Figure 3. 

Figure 3 shows that Rm/Rout ratios > 0.5 gen- 
erally reproduce the results, depending on the SFE 
and Rout values adopted. The very high SFE val- 
ues, such as 50%, always allows a large fraction of 
the ultraviolet photons out of the envelope because 
of its low optical depth. The dumpiness will also 
have an effect on the fraction of UV light leakage, 
and generally less light will leak from envelopes 
with smaller clumpy fractions. 

Since the radio sources studied in Reines et 
al. (2008) and Johnson et al. (2009) are obscure, 
however, their findings are probably for low-SFE, 
moderately evolved regions where the optical ex- 
tinction for a corresponding model is about Ay ~ 
30 — 50 on average. The UV light leakage would 
therefore take place along those sightlines exhibit- 
ing minimal extinction. Ay < 1. 
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Fig. 3. — The fraction of UV light that escapes from 
the 99% clumpy embedding envelope as a function of 
Rijj/Rout- The horizontal line is at a value of 0.5, sig- 
nifying 50% of the UV flux is not being reprocessed by the 
dust, and is only meant to help guide the eye. The solid 
line is for a SFE of 5%, dotted line for 10%, dashed line for 
15%, dash-dot line for 25%, and dash-dot-dot for 50%. 



3.3. Infrared Variation with Viewing An- 
gle 

The clumpy dust distribution creates variations 
in infrared luminosity and certain spectral fea- 
tures depending on the viewing angle. In a gen- 
eral sense, clumpier media show the most varia- 
tion in spectral features and luminosity compared 
to smoother media. 

Figure 4 demonstrates how the infrared lumi- 
nosity derived from a random sightline will be in- 
correct due to a clumpy envelope. The infrared 
luminosity inferred from a single sightline obser- 
vation, plotted as stars in the figure and defined 
as the 3-1000/zm integrated luminosity, can vary 
from roughly half to nearly twice the true in- 
frared luminosity for a clumpy dust distribution. 
The true infrared luminosities are plotted as dia- 
monds and are the infrared luminosities summed 
over all viewing angles for a model. For smooth 
dust distributions with thick envelopes, the input 
central source luminosity is equal to the true in- 
frared luminosity. The infrared luminosities are 
compared to the stellar cluster's input luminosity 
of 1.604 X IO^Lq, plotted as the solid horizontal 
line. 

If a clump is along the line of sight in front 
of the central source, the optical depth will ap- 
pear to be high and the 9.8 and 17^m silicate 
absorption features will be deeper in the infrared 



spectrum. However, clumps that are behind the 
central source will reflect light, and silicate emis- 
sion features will therefore be evident. This is il- 
lustrated in Figure 5, which shows that the fea- 
tures of the infrared spectrum of a source inside 
a clumpy envelope are highly sightline-dependent. 
The average spectrum is also shown for compari- 
son. 
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Fig. 4. — The computed infrared luminosity (3-lOOO^m) 
plotted versus the fraction of dust that is in clumps for an 
optically thick (top) and an optically thin (bottom) model. 
Kout is 25pc for both panels. The horizontal line is the 
input stellar cluster luminosity, which matches the Lju 
for optically thick models. The Lm computed along eax:h 
sightline, assuming a spherical, smooth distribution, is plot- 
ted with stars, while the Lju values determined by adding 
up all sightlines for a model are shown as diamonds (this 
is the correct infrared luminosity). The clumpier models, 
as well as showing a larger dispersion in computed infrared 
luminosity values, also exhibit lower infrared luminosities 
in general, because the clumpy envelope allows ultraviolet 
photons to exit the envelope (see Figure 3). 



3.4. Spectral Energy Distribution Proper- 
ties 

The output spectral energy distributions all 
generally have four components: (1) an extincted 
stellar spectrum from the central source; (2) ther- 
mal dust emission that peaks between 40 and 
70/Ltm; (3) silicate emission or absorption features; 
and (4) the PAH features visible in the mid-IR. 
Figure 6 shows the smooth dust SEDs for a geo- 
metric sequence as an example of how these four 
main components are expected to change as the 
cluster evolves. The geometric sequences for all 
SFE and Rout values are given in Figures 7 and 
8 for the smooth and 99% climipy dust distribu- 
tions. In Figure 8, the gray spectra are individual 
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Fig. 5. — The variation of the SED for a given model 
due to the viewing angle of the observer. For this example, 
the model shown has Rout = 25 pc, Rin = 3 pc, a clumpy 
fraction of 0.90, and a 10% SFE. The mean SED is shown 
with a black line, and the near-IR and Spitzer observing 
bands are indicated with diamonds, with AA for each band 
given at the bottom. 

sightline SEDs for each Ri„ value, while the colors 
are the sightline-averaged values^. 

For the most dense models, including low SFE 
values and small values, the starlight is largely 
absorbed by the envelope. For clumpy envelopes, 
many sightlines allow optical and UV photons to 
pass through so that a stellar continuum is visible, 
while the clumps have very high extinction val- 
ues. For smooth dust distributions the starlight 
is largely absorbed along all sightlines. The Ay 
values for smooth distributions compared to the 
average values and ranges found in clumpy solu- 
tions, are listed in Table 2 for the Ri„ = O.lpc 
models. The average Ay values for evolved stages 
in the geometric sequence are relatively small for 
all models (roughly less than 2), hence the spectra 
of clumpy and smooth models at high SFE value 
and large inner radii are similar. 

The wavelength of the far-IR peak of the ther- 
mal grain emission depends on the predominant 
temperature of the dust in a model envelope. For 
the high column density envelopes, which have low 
SFE values and small R™ values, the UV light 
is absorbed in the inner portions of the envelope 
so that much of the envelope remains quite cold. 
Therefore, low SFE values and small R^ values 
tend to produce infrared peaks at about 60fim. 
For the optically thin models with high SFE val- 
ues and large inner radii, a large proportion of the 



dust is being heated by the cluster, so the far-IR 
peak moves to shorter wavelengths. 

Silicate absorption is directly proportional to 
the optical depth, so that the deepest silicate fea- 
tures appear in the largest column-depth models. 
Silicates in emission, as was described in § 3.3, 
appear in clumpy models where sightlines include 
dust clumps illuminated behind the central source, 
which re-emit absorbed starlight as 10/im light 
that is unblocked by any intervening clouds. 

The PAH features are almost ubiquitous in the 
model SEDs. The early stages in the geometric 
sequence appear to have reduced emission from 
the PAHs for two reasons. The first reason is 
that for the high density models, there is very lit- 
tle volume from which the PAHs can be excited; 
PAHs are only excited by UV photons to a depth 
of Ay = 1.5. The second reason is because the 
thermal hot dust component dominates the mid- 
IR flux for the early stage models, which lowers the 
measured equivalent widths of these features. For 
later stages on the geometric sequence, the PAH 
emission is relatively constant along the sequence 
because as the ultraviolet flux on the inner surface 
of the envelope decreases as r~^, the surface area 
available to excite PAH molecules increases as r^. 




wQvelergth (urn) 



Fig. 6. — The spectra shown are along the geometric 
sequence in Table 1 for Rout = 50pc and a smooth dust 
distribution. As Ri„ increases towards Rout, the silicate 
features, most noticeably the one at ~ 9.8/im, become shal- 
lower. The PAH features are about the same throughout 
the sequence. The envelopes with small inner radii are more 
opaque to short-wave radiation and reprocess it to longer 
wavelengths. Therefore less flux between 1-5 /im is mea- 
sured and more flux at 100 ^m, shifting the peak of the 
SED to the red. The stellar light is less extincted at larger 
Hin values. 



The suite of model SEDs and photometric bands along all 
viewing angles are available via D. G. Whelan's webpage: 
http : / /www. astro.virginia.edu/ ~ dwwlv / S D jmodels / 
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Fig. 7. — The infrared SEDs of all of the geometric se- 
quences using a smooth dust distribution are shown; the 
range in Ri„ values from small to large are as in Figure 6. 



Fig. 8. — The infrared SEDs of all the geometric se- 
quences using a 99% clumpy dust distribution, with the 
sightline-average data shown in colors as in Figures 6 and 
7 and the individual sightlinc data shown in gray scale. 



3.5. Infrared Colors as Tracers of the In- 
put Parameters 

The model SEDs have been convolved with nu- 
merous filter sets for comparison to photometric 
observations. These include the four IRAS bands 
(12, 25, 60 and 100/xm), three standard near-IR 
bands (J, H, and K), Spitzer IRAC bands (3.6, 
4.5, 5.8, and 8.0/xm), Spitzer MIPS bands (24, 70, 
and 160/Ltm), Herschel PACS bands(70, 100, and 
160//m) and 77er.se/ie/SPIRE bands (250, 350, and 
500//m). 

As an example of how colors can be used to plot 
the data, we show the IRAC colors for the Rout = 
25pc and SFE= 10% models for all dumpiness 
values and all Rm values in Figure 9. There are 
degeneracies in the model average values at late 
evolutionary stages, and the individual sightlines 
make the interpretation of data extremely difficult. 



Because the clumpy dust distributions create 
degeneracies in the SEDs, diagnostics to recover 
the input parameters were developed using the 
smooth dust distribution models first. In Fig- 
ure 10 the [3.6]- [24] color is plotted against the 
Rm/Rout- Red, green, and blue are used to desig- 
nate models with different Rout values (red, green, 
and blue for 100, 50, and 25pc, respectively), and 
the plotting symbols are for the five different SFE 
values considered. Although there is significant 
overlap between the colors, the tracks are separate 
over most of the Rin/Ro«t values for each individ- 
ual SFE value. The [70]-[160] color can be used to 
roughly determine the SFE value as shown in Fig- 
ure 11; the bars represent the model dispersions. 
There is significant overlap between dispersion in 
SFE values of 10% to 25%, but [70]- [160] values 
of > 0.9 represent very low SFE values (< 10%), 
values between ^ —1.0 and —1.5 represent SFE 
values between 10% and 25%, and values > — 1.6 
will be for SFE of 50% or greater. 

Colors that can best be used to recover the 
input parameters for the entire, clumpy suite of 
models were also investigated. Given that the 
value of the input parameters and the viewing 
angle can both significantly change an embedded 
source's SED and colors, there is a need to deter- 
mine which, if any, photometric measurements can 
be used to reliably constrain the physical geome- 
try. For a given input parameter (SFE, dumpi- 
ness, Rm, or Rout), the following calculation was 
made: for each color, the mean color and standard 
deviation was measured with the input parameter 
fixed and everything else (i.e. the other input pa- 
rameters and all viewing angles) variable. This 
was done at each fixed input parameter value (e.g. 
for each of the five different SFE values). Finally 
the difference in the means is divided by the great- 
est of the standard deviations - this is a measure 
of how much the input parameter affects the color, 
compared to the other input parameters and the 
viewing angle ambiguities. 

This analysis is related to a Principal Compo- 
nent Analysis, but adjusted to our goal of finding 
the minimal set of colors with maximal physical 
diagnostic power. Rather than solving freely for 
eigenvectors of the model set in color space, which 
may not directly correspond to the physical vari- 
ables, we hypothesize that a principal component 
exists for which the eigenvalue would directly cor- 
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respond to a physical parameter, and then mea- 
sure the projection of that hypothetical compo- 
nent onto each of our color axes. This process 
allows of course that there is no such 'physically 
diagnostic principal component,' in which case the 
variation of all colors with the physical parameter 
would be small compared to the standard devia- 
tion in the color due to other causes. 

Figures 12 to 14 are plots of infrared color ver- 
sus infrared color where the grayscale shows the 
value of the quotient for determining the best diag- 
nostic color described above. Values in the figures 
near zero indicate little spread in the average color 
for each input parameter value, while values at the 
maximum have relatively large spreads in average 
color values across the input parameter values, as 
well as small standard deviations from the mean. 
These are therefore the best colors for separat- 
ing input parameter values: (1) [70]- [160] for SFE 
(Figure 12); (2) [3.6]-[5.8] for fraction of dust in a 
clumpy distribution (Figure 13); and (3) [8.0]- [24] 
for both Rout and the Rin/Rout ratio (Figure 14). 

Figures 15, 16, and 17 plot these most useful 
colors' average values and standard deviations for 
the whole suite of models versus the input param- 
eters. For SFE (Figure 15), as with the smooth 
dust distribution models, there is overlap in the 
[70]- [160] color between the various SFE values. 
The overlap is still most significant between 10 and 
25%, for which color values between —1.5 and —0.9 
would all be valid. Color values above —0.9 match 
the 5% SFE, and color values below —1.5 match 
the 50% SFE. For the smooth dust fraction, within 
the error bars in the [3.6]-[5.8] versus smooth dust 
fraction plot (Figure 16) the color values are still 
degenerate, and using thermal radio data (such as 
was presented in Reines et al. 2008) is a more 
reliable metric for determining the fraction of ul- 
traviolet light lost from an embedded star cluster. 
However, because the standard deviations in color 
of the clumpiest fractions are small compared to 
the smooth models (due to the large differences 
in the input parameter values being explored), a 
large, statistical study of embedded star clusters, 
for which average color values and standard de- 
viations could be computed, could potentially be 
used to discern between the different smooth dust 
fraction values. For Rm/Rout (Figure 17), colors 
are still intractable for small Rm/Rout values. For 
'Rin/Rout > 0.4, the [8.0]-[24] can be used both to 



distinguish between different Rin/Rout values and 
to distinguish between the different Rout values. 



SFE = 10% 




2 f . ^ ; 

2 4 6 8 10 

[3.6]-[5,8] 



Fig. 9. — Using the mid-IR Spitzer/IRAC colors, the 
geometric sequence for Rout = 25pc and SFE= 10% is 
shown for smooth models (black line), 10% clumpy models 
(purple), 50% clumpy models (blue), 90% clumpy models 
(green), and 99% clumpy models (red). The geometric se- 
quences begin towards the bottom with Ri„ = O.lpc and 
follow the arrows around the turn-around to their termina- 
tion at Ri„ = 24pc. The solid lines and diamonds trace out 
the sightline-averaged colors, and the dots represent indi- 
vidual sightlines. While even the individual average tracks 
often overlap, the different sightlines add a significant un- 
certainty to any interpretation of these data. 
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Fig. 10. — The [3.6]-[24] color is plotted versus the ratio 
of the inner and outer envelope radii for the three different 
outer radius values considered, 25, 50, and lOOpc (shown in 
blue, green, and red, respectively) using the smooth dust 
distribution models only. If the star formation efficiency 
(SFE) can be determined using the [70]-[160] color as shown 
in Figure 11, then the Ri„/Rout value can be estimated 
using this plot. 




[3.6][4.5][5.8][8.0][24] [70][ 1 60l250l500] 
band 



Fig. 12. — Determining the best color (y-axis magnitude 
minus x-axis magnitude) for recovering the star formation 
efficiency (SFE). In order to determine the best color for 
discriminating between different SFE values, this calcula- 
tion was made: the spread in the mean color values be- 
tween the different SFE values was divided by the greatest 
standard deviation of each SFE value was computed. The 
maximum resulting value is the most discriminating color 
(see §3.5). Each color is the y-axis magnitude minus x-axis 
magnitude, and the best color was found to be [70]-[160]. 
How well this 'best' color helps to determine the SFE is 
shown in Figure 15, and discussed fully in §3.5. 




Fig. 11. — The [70]- [160] color is used to separate differ- 
ent SFE values for the models with a smooth dust distribu- 
tion. The bars shown are the dispersion and there is a fair 
bit of overlap between the models. One can only expect, 
therefore, to be able to roughly determine the SFE for the 
smooth models. 



Fig. 13. — Determining the best color for recovering the 
fraction of dust in a smooth distribution in the envelope. 
Using the criterion described in §3.5, the best color for 
recovering the percentage of dust that is smooth was found 
to be [3.6]-[5.8]. However, as shown in Figure 16, even this 
'best' color cannot be used to discriminate between the 
average values. 
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[3.6][4.5][5.8][8.0][24] [70][ 1 60l250l500] 
bond 



Fig. 14. — Determining the best color for recovering the 
Rout ratio. The best discriminator found for Rout, which 
also well separates and Rin/Rout ratios, is [8.0]-[24], and is 
shown in Figure 17. 
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Fig. 15. — The [70]-[160] color that can best be used 
to recover the SFE value. The diamonds are the average 
values over the entire dataset, and the error bars represent 
the one-cr errors at each SFE value. 
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Fig. 17. — For envelopes with Ri„/Rout > 0.4 the [8.0]- 
[24] color can be used to determine both the Rout value, 
because they are well separated from each other, and the 
Rin/Rout ratio. 
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Fig. 16. — The [3.6]-[5.8] color that was found to best 
recover the fraction of dust that is smoothly distributed in 
the envelope. Given the l-cr error bars for the smoothest 
models (i.e. there is such a great diversity in this color 
due to the input parameters) it is only possible to sepa- 
rate the smooth from the clumpy models in a statistical 
sample where the standard deviation in the color could be 
measured. 
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4. COMPARISONS TO OBSERVATIONS 

In the Milky Way, individual embedded massive 
stars are surrounded by what have been identified 
as "Ultracompact Hii Regions" (UCHiis) , and it is 
hypothesized that the constituent massive stars in 
embedded SSCs will also be surrounded by UCHii 
regions. For this reason we compare our models 
to a sample of resolved embedded massive stars 
in order to test whether our models have repro- 
duced, to first order, the necessary features of em- 
bedded massive stars. The sample is an IRAS 
sample studied in Wood & Churchwell (1989), 
in which it was found that 60% of the brightest 
IRAS sources (> 10^ Jy at lOO/im) in the color 
range of their galactic survey are UCHii regions. 
Wood & Churchwell (1989) also found that galac- 
tic UCHii regions strictly obey a set of color crite- 
ria in the infrared with log{FQQ^m/ Fi2fj,m) > 1-30 
and log{F25^m/ Fi2^m) > 0.57, while very few 
other types of objects had IRAS colors fitting these 
criteria. Therefore, these color criteria appear to 
be relatively robust for identifying UCHii regions, 
and we might expect natal SSCs to obey these cri- 
teria as well. 

We have convolved our model results with the 
IRAS filter profiles, and Figure 18 shows the 
IRAS colors from a subset of our models com- 
pared to field objects (plus signs) and UCHii re- 
gions (diamonds) from the Kurtz, Churchwell, & 
Wood (1994) survey (which includes the Wood 
& Churchwell (1989) sample). In this figure, we 
focus on only the family of models with a 10% 
stellar mass component and a relative fraction of 
clumpy dust of 90%, although the other model 
families have IRAS colors that are included in the 
range represented by the models shown. Mod- 
els with RoMt — 25 pc are shown as red points. 
Rout = 50 pc as green points, and Rout = 100 pc 
as purple points. The values of Ri„ used are those 
shown in Table 1, and all sight-lines are shown 
for each set of model parameters, which produces 
additional scatter. 

There is good agreement between our models 
and the UCHiis sample shown in Figure 18. There 
are a few of points to lower left and upper right 
that are not traced by our models. The lower left 
points have very hot dust in the inner portions 
of their envelopes, and might be modeled with a 



smaller Kin value that puts the hot dust closer to 
the heating source. The points to upper right seem 
to be an extension of the basic trend of increasing 
Rin- It is therefore possible that we could repro- 
duce these UCHiis if we simply assumed a lower 
mass cutoff to our cluster which would lower the 
UV flux incident on the inner envelope. 




1 ID 100 

Fig. 18. — The IRAS colors from a subset of our models 
assuming 10% SFE and clumpy fraction = 0.9 compared 
to field objects (plus signs) and UCHll regions (diamonds) 
from the Kurtz, Churchwell, & Wood (1994) survey. Mod- 
els with Rout = 25 pc are shown as red points, Rout = 50 pc 
as green points, and Rout = 100 pc as purple points. The 
values of Rj^ used are those shown in Table 1, and all sight 
lines are shown for each set of model parameters. 

5. KNOWN LIMITATIONS 

Here we discuss some of the limitations of the 
models as a guide for future improvements. 

(1) The models assume that the central star 
cluster is a point source. While optical mea- 
surements have discovered half-light radii of super 
star clusters of 1.7pc on average (de Marchi et al. 
1997), which is significantly smaller than the as- 
sumed outer envelope radii, it is nevertheless more 
realistic to place several ultracompact Hii regions 
at the center of the envelope (see Figure 1). The 
intracluster dust would be heated at early evo- 
lutionary stages if multiple, smaller clusters were 
used, creating a hotter dust component than is 
seen with the current models. 

(2) Variations in metallicity have not been con- 
sidered. For small changes in metallicity (i.e. 
down to half or third solar metallicity), we ex- 
pect the dust-to-gas ratio to scale with metallicity. 
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However, more extreme environments, like blue 
compact dwarf galaxies, where metallicities are of- 
ten much lower than this, show dust in excess to 
what is expected from their metallicity alone (T. 
X. Thuan, private communication). Accounting 
for very low metallicities is therefore a complicated 
issue that is probably best handled on an individ- 
ual basis. 

(3) The evolutionary sequences presented in 
this paper do not allow the central star cluster 
to evolve. As described in §2.5, the time between 
cluster formation and the dispersal of the embed- 
ding envelope is short, on the order of 3 — 4 Myr. 
However, the dispersal rate is not well known, 
so we have not made any attempt to place time 
stamps on the evolutionary sequence presented 
here. Hydrodynamical models are needed to an- 
swer this problem, and several hydrodynamical 
and semi-analytical studies of super star clusters 
already exist (e.g. Tenorio-Tagle et al. 2005, 2007; 
Wiinsch et al. 2007, 2008). The rate at which 
the envelope is dispersed could potentially be an- 
swered by similar studies of a super star cluster's 
embedded phase. 

(4) This work has concentrated solely on the 
dust emission and absorption in embedded super 
star clusters. Inclusion of nebular line emission 
in the model spectra could be useful if line diag- 
nostics could be identified that break the model 
degeneracies. 

(5) The models do not include an interstellar 
radiation field (ISRF) incident on the outside of 
the dusty envelope. An ISRF would heat the dust 
on the outside of the envelope and provide an ad- 
ditional component to the starlight visible in the 
UV/optical/near-IR regime. This could affect the 
[3.6]- [5.8] color used to distinguish between dif- 
ferent fractions of chimpy dust. The heating of 
the outer envelope could affect the far-IR fluxes 
as well. 

6. SUMMARY 

We present SEDs and colors of embedded SSCs 
created using spherical three-dimensional models. 
By varying the input parameters according to a 
series of evolutionary sequences, we have created 
a suite of models that can be used to constrain 
the evolutionary state of an embedded super star 
cluster. The main conclusions of the study are: 



(1) A hierarchically clumped rncidium is suitable 
for recreating the porous environments observed 
around embedded super star clusters (e.g. Johnson 
et al. 2009); 

(2) The infrared luminosity derived from a sin- 
gle sightline observation of a clumpy envelope can 
be wrong by as much as a factor of two from the 
true value; 

(3) The infrared SED also varies with sight- 
line in these clumpy models, which has an impact 
on the near- and mid-infrared colors and magni- 
tudes, the strength of the observed silicate fea- 
tures, and the dust continuum measured at these 
wavelengths; 

(4) For the smooth dust distribution models, 
the evolutionary sequences that begin with mid- 
infrared obscure envelopes (Ay ^ 50) are marked 
by a gradual decline in the silicate absorption fea- 
tures at 9.8 and 17/zm and a corresponding in- 
crease in the visual and ultraviolet flux as the clus- 
ter envelope evolves. Those sequences that begin 
with infrared-transparent envelopes {Ay < 50) in- 
stead have a predominant hot dust component and 
silicates in emission at early stages that eventually 
both fade away as the inner envelope radius moves 
outward. The clumpy envelope acts to confuse 
these general trends, making it harder to deter- 
mine envelope properties. 

(5) Several diagnostic colors were found to con- 
strain the envelope properties. The Spitzer MIPS 
[70]- [160] color is found to be a good diagnostic of 
the star formation efficiency, particularly at sepa- 
rating very low and high values (such as 5% and 
50%) from more moderate values (between 10% 
and 25%). The [3.6]-[5.8] color can be used to 
determine the fraction of clumpy dust in the en- 
velope for large samples of embedded super star 
clusters, but not for individual sources (see §3.5 
for details). In order to determine the Ri„ and 
Rout values, the [8.0]-[24] color can be used for 
Rjn/Rout ^ 0.4. Below this value the data is de- 
generate for all colors. 

(6) The model IRAS colors trace the same area 
of color space as ultracompact Hii regions, the 
Galactic analogues to extragalactic, embedded sni- 
per star clusters, suggesting that the models will 
also be useful when data of resolved, embedded 
super star clusters become available. 
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Table 1 
The geometric sequences 



Rout = 25pc 


R-out = 50pc 


Rout = lOOpc 


Rin 


Rin 


Rin 


O.lpc 


O.lpc 


O.lpc 


0.5pc 


0.5pc 


l.Opc 


l.Opc 


l.Opc 


5.0pc 


2.0pc 


5.0pc 


lOpc 


3.0pc 


lOpc 


15pc 


e.Opc 


15pc 


25pc 


g.opc 


20pc 


35pc 


12pc 


25pc 


45pc 


15pc 


30pc 


55pc 


18pc 


35pc 


65pc 


21pc 


40pc 


75pc 


24pc 


45pc 


95pc 



Note. — Each sequence is run using clumpy dust fractions of 0.0, 0.1, 0.5, 0.9, and 0.99, and SFEs of 
5%, 10%, 15%, 25%, and 50%, resulting in a total of 900 models. Stellar cluster mass of W^Mq, source 
luminosity ~ 1.6 x IO^Lq, Salpeter IMF 1 - lOOM©, age 1 Myr, Zq. 



Table 2 

A table of average Ay VALUES AND RANGES FOR THE Rj„ = O.lPC MODELS 



SFE(%) 




Smooth 


10% clumpy 


50% clumpy 


90% clumpy 


99% clumpy 


5 


25 


490 


474 (441-569) 


411 (245-885) 


348 (49.1-1200) 


334 (4.92-1270) 


5 


50 


123 


119 (111-143) 


103 (61.4-221) 


87.0 (12.3-300) 


83.4 (1.24-318) 


5 


100 


30.8 


29.8 (27.6-35.7) 


25.8 (15.4-55.4) 


21.8 (3.09-75.1) 


20.9 (0.310-79.5) 


10 


25 


232 


225 (209-270) 


195 (116-419) 


165 (23.2-569) 


158 (2.34-603) 


10 


50 


58.2 


56.3 (52.4-67.5) 


48.8 (29.1-105) 


41.2 (5.83-142) 


39.5 (0.587-151) 


10 


100 


14.6 


14.1 (13.1-16.9) 


12.2 (7.30-26.2) 


10.3 (1.47-35.6) 


9.88 (0.147-37.7) 


15 


25 


146 


142 (132-170) 


123 (73.2-264) 


104 (14.6-358) 


99.5 (1.48-379) 


15 


50 


36.7 


35.5 (33.0-42.5) 


30.7 (18.3-66.1) 


26.0 (3.68-89.6) 


24.9 (0.369-94.8) 


15 


100 


9.18 


8.88 (8.27-10.7) 


7.69 (4.60-16.5) 


6.50 (0.925-22.4) 


6.23 (0.0925-23.7) 


25 


25 


77.5 


75.0 (69.7-89.9) 


64.9 (38.7-140) 


54.9 (7.76-190) 


52.7 (0.781-201) 


25 


50 


19.4 


18.8 (17.5-22.5) 


16.3 (9.71-35.0) 


13.7 (1.95-47.4) 


13.2 (0.196-50.2) 


25 


100 


4.87 


4.71 (4.38-5.65) 


4.08 (2.44-8.76) 


3.44 (0.490-11.9) 


3.30 (0.0490-12.6) 


50 


25 


25.8 


25.0 (23.2-30.0) 


21.7 (12.9-46.6) 


18.3 (2.59-63.2) 


17.6 (0.260-67.0) 


50 


50 


6.48 


6.27 (5.83-7.52) 


5.43 (3.25-11.7) 


4.59 (0.652-15.8) 


4.40 (0.0652-16.7) 


50 


100 


1.63 


1.58 (1.47-1.89) 


1.36 (0.816-2.93) 


1.15 (0.163-3.96) 


1.10 (0.0163-4.20) 



Note. — The Ri„ value is O.lpc for all of these values. The range in values depending on sightline are in parentheses for 
the clumpy models, with the average values shown first. 
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